I want to determine if there is any change in the peak-gene linkage for differentially expressed genes across conditions. While the target genes are going to be primarily from the cytotoxic data, I do want to check all the populations to see if its a more common thing.
suppressPackageStartupMessages(library(ggplot2))
suppressPackageStartupMessages(library(dplyr))
suppressPackageStartupMessages(library(Seurat))
suppressPackageStartupMessages(library(Signac))
results<-readRDS("~/gibbs/DOGMAMORPH/Ranalysis/Objects/20230628Linkpeakobj.rds")
#specifically made this for my morphine stuff, if someone else uses it in the future they should pay attention to the number of plots and the order. There are ways to make it more general but I opted to make it as close to what I need as possible
PlotLinked<-function(objlist, gene,idents.plot=c(""), up=10000, down=10000, expression.assay="RNA", assay="ATAC" ){
if(idents.plot[[1]]==""){idents.plot<-names(table(Idents(objlist[[1]])))}
p1<-CoveragePlot(objlist[[3]], region=gene, features=gene, extend.downstream = down, extend.upstream = up, assay = assay,expression.assay = expression.assay, idents = idents.plot)
p2<-CoveragePlot(objlist[[4]], region=gene, features=gene, extend.downstream = down, extend.upstream = up, assay = assay,expression.assay = expression.assay, idents = idents.plot)
p3<-CoveragePlot(objlist[[1]], region=gene, features=gene, extend.downstream = down, extend.upstream = up, assay = assay,expression.assay = expression.assay, idents = idents.plot)
p4<-CoveragePlot(objlist[[2]], region=gene, features=gene, extend.downstream = down, extend.upstream = up, assay = assay,expression.assay = expression.assay, idents = idents.plot)
p5<-CoveragePlot(objlist[[5]], region=gene, features=gene, extend.downstream = down, extend.upstream = up, assay = assay,expression.assay = expression.assay, idents = idents.plot)
p6<-CoveragePlot(objlist[[6]], region=gene, features=gene, extend.downstream = down, extend.upstream = up, assay = assay,expression.assay = expression.assay, idents = idents.plot)
return((p1+p2)/(p3+p4)/(p5+p6)+ patchwork::plot_layout(nrow = 3, byrow = FALSE) )
}
PlotLinked(results, gene = "CD4")
## Warning: Removed 85 rows containing missing values (`geom_segment()`).
## Removed 85 rows containing missing values (`geom_segment()`).
## Removed 85 rows containing missing values (`geom_segment()`).
## Removed 85 rows containing missing values (`geom_segment()`).
## Removed 85 rows containing missing values (`geom_segment()`).
## Removed 85 rows containing missing values (`geom_segment()`).
PlotLinked(results,"OPRM1", up = 1e5, down = 1e5 )
## Warning: Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
PlotLinked(results,"OGFR", up = 1e5, down = 1e5 )
## Warning: Removed 13 rows containing missing values (`geom_segment()`).
## Warning: Removed 13 rows containing missing values (`geom_segment()`).
## Removed 13 rows containing missing values (`geom_segment()`).
## Removed 13 rows containing missing values (`geom_segment()`).
## Removed 13 rows containing missing values (`geom_segment()`).
## Removed 13 rows containing missing values (`geom_segment()`).
PlotLinked(results,"IFI44L", up = 1e5, down = 1e5, idents.plot = c("Effector_Th1_CD4_and_CD8","CD8_TEMRA") )
## Warning: Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
PlotLinked(results,"PLCG2", up = 1e5, down = 1e5 , idents.plot = c("Effector_Th1_CD4_and_CD8","CD8_TEMRA") )
## Warning: Removed 29 rows containing missing values (`geom_segment()`).
## Warning: Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
PlotLinked(results,"BRD1", up = 1e5, down = 1e5 , idents.plot = c("Effector_Th1_CD4_and_CD8","CD8_TEMRA") )
## Warning: Removed 5 rows containing missing values (`geom_segment()`).
## Warning: Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
PlotLinked(results,"AKT3", up = 1e5, down = 1e5 , idents.plot = c("Effector_Th1_CD4_and_CD8","CD8_TEMRA") )
## Warning: Removed 60 rows containing missing values (`geom_segment()`).
## Warning: Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
PlotLinked(results,"CCL5", up = 1e5, down = 1e5 , idents.plot = c("Effector_Th1_CD4_and_CD8","CD8_TEMRA") )
## Warning: Removed 22 rows containing missing values (`geom_segment()`).
## Warning: Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
PlotLinked(results,"SPOCK2", up = 1e5, down = 1e5 , idents.plot = c("Effector_Th1_CD4_and_CD8","CD8_TEMRA") )
## Warning: Removed 72 rows containing missing values (`geom_segment()`).
## Warning: Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).
PlotLinked(results,"NKG7", up = 1e5, down = 1e5 , idents.plot = c("Effector_Th1_CD4_and_CD8","CD8_TEMRA") )
## Warning: Removed 3 rows containing missing values (`geom_segment()`).
## Warning: Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
## Removed 3 rows containing missing values (`geom_segment()`).
PlotLinked(results,"IFI44L", up = 1e5, down = 1e5)
## Warning: Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
## Removed 11 rows containing missing values (`geom_segment()`).
PlotLinked(results,"PLCG2", up = 1e5, down = 1e5)
## Warning: Removed 29 rows containing missing values (`geom_segment()`).
## Warning: Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
## Removed 29 rows containing missing values (`geom_segment()`).
PlotLinked(results,"BRD1", up = 1e5, down = 1e5)
## Warning: Removed 5 rows containing missing values (`geom_segment()`).
## Warning: Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
## Removed 5 rows containing missing values (`geom_segment()`).
PlotLinked(results,"AKT3", up = 1e5, down = 1e5)
## Warning: Removed 60 rows containing missing values (`geom_segment()`).
## Warning: Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
## Removed 60 rows containing missing values (`geom_segment()`).
PlotLinked(results,"CCL5", up = 1e5, down = 1e5)
## Warning: Removed 22 rows containing missing values (`geom_segment()`).
## Warning: Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
## Removed 22 rows containing missing values (`geom_segment()`).
PlotLinked(results,"SPOCK2", up = 1e5, down = 1e5)
## Warning: Removed 72 rows containing missing values (`geom_segment()`).
## Warning: Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).
## Removed 72 rows containing missing values (`geom_segment()`).